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Abstract 

For a system of ./V qubits, spanning a Hilbert space of dimension d = 2 , it is known that 
there exists d + 1 mutually unbiased bases. Different construction algorithms exist, and it is 
remarkable that different methods lead to sets of bases with different properties as far as separability 
is concerned. Here we derive the four sets of nine bases for three qubits, and show how they are 
unitarily related. We also briefly discuss the four-qubit case, give the entanglement structure of 
sixteen sets of bases, and show some of them, and their interrelations, as examples. The extension 
of the method to the general case of N qubits is outlined. 

PACS numbers: 3.65.Ta, 03.65.Ud, 03.65. Ca 
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I. INTRODUCTION 



Every quantum system is associated with some state (pure or mixed) in a Hilbert space. It 
is possible to ascertain this quantum state by performing a series of measurements on an en- 
semble consisting of many identically members. Each measurement will modify the measured 
ensemble member in such way that it, in general, not possible to get any additional informa- 
tion about the original state. Several techniques, such as state tomography , maximum 
likelihood jjj U, IJ , and maximum-entropy methods |(| (or combinations thereof j?j|), have 
been devised for state estimation. 

When the Hilbert space is finite, it has been shown that the optimal approach to get the 
information is related to a special set of states that are "mutually unbiased" J^, 0, , for 
which the uncertainty spread of the inferred state is minimized. Note, however, that we are 
ignoring more general measurements, such as joint measurements on all the members of the 



more gener; 



l4 15^ , which surpass the ability of a priori 



ensemble IllL I12L 1131 1 or adaptive measurements 
fixed, single-copy measurements. 

Let us denote basis sets by A = 1,2, . . . and states within a basis by \A, a), with a = 
1,2, ... ,d, d being the dimension of the Hilbert space. We recall that two bases \A, a) and 
\B,b) are said to be mutually unbiased bases (MUBs) if a system prepared in any element 
of A has a uniform probability distribution of being found in any element of B, that is 



\{A,a\B,b)\ 2 = l/d, 



where orthonormality among states of the same basis is assumed. These MUBs are central 



to the formulation of the discrete Wigner function 



been used in cryptographic protocols 
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. They have also have also 
2l| . due to the complete uncertainty about the 
outcome of a measurement in some basis after the preparation of the system in another, if the 
bases are mutually unbiased. MUBs are also used for quantum error correction codes [22], |2J| 



and recently they have also found uses in 



solution to the "mean king problem" 



quantum 
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27. 



ame theory, in particular to provide a 



It has been shown that the maximum number of MUBs can be at most d+ 1 [8j]. Actually, 
it is known that if d is prime or power of prime, the maximal number of MUBs can be 
achieved . Remarkably though, there is no known answer for any other values of d, 

although there are some approaches to try to find the solution to this problem in some simple 



2 



cases, such as d = 6 or when d is an nonprime integer squared |30|, |3jJ, |32|]. Recent works 
have suggested that the answer to this question may well be related with the non-existence of 
finite projective planes of certain orders (33, 34] or with the problem of mutually orthogonal 



Latin squares in combinatorics 
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Experimental quantum information and computation have already moved from single- 
qubit protocols to several qubits (at present around eight 3?j]), so there is a need to extend 
our knowledge, especially about entanglement properties of several qubits. This also includes 
extensions of measurement techniques of systems with more than two qubits. Therefore, a 
new problem related with MUBs naturally appears, namely that for more than two qubits, 
different MUB structures exist, where the word "structure" refers to the entanglement prop- 
erties of the bases. We are already aware of the existence of three MUB structures for three 
qubits nll^. In this paper we will show that, in fact, there exists exactly four different 
MUB structures in this space. We will also show how they are interrelated. For the experi- 
mentalist, this information is very important, because the complexity of an implementation 
of two or more MUBs will, of course, greatly depend on how many of the qubits that need 
to be entangled. We will also briefly discuss the four-qubit case and show that in this space 
there exist sixteen different MUB structures. We will exhibit which they are and derive 
some of them explicitly. It is then possible to continue and analyze the general ./V-qubit 
MUBs much in the same manner, although, for brevity and simplicity, we will stop at four 
qubits. 



II. MUBS FOR ONE AND TWO QUBITS 

Because states belonging to the same basis are usually taken to be orthonormal, to study 
the property of "mutually unbiasedness" it is possible to use either mutually unbiased bases 
or the operators which have the basis states as eigenvectors. We thus need d 2 — 1 operators 
to obtain the whole set of states. In the case of prime and power of a prime dimensions this 
set of operators can be collected as d + 1 sets of d — 1 commuting operators, which is related 
with the grading of a Lie algebra |39:]. In the case of iV-qubits (2^- dimensional case) we 
need iV commuting operators to define uniquely a pure state Q,Q- 

In finite dimensional systems is also possible to define a discrete phase space, and when 
the dimension of the system is either a prime or a power of prime the phase space is a 
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finite geometry. The above operators are related with translations in this phase space and 
(without a phase factor) they are the so-called displacement operators, which satisfy the 
covariant property of the discrete Wigner function defined there 18]. 

The two-dimensional Hilbert vector space (one qubit) is spanned, e.g., by the two or- 
thonormal eigenvectors of the spin 1/2 observable a z , which will be used, in the following, 
as our computational basis. 

In this Hilbert space, the MUB set of 2 2 — 1 = 3 bases is given by the eigenvectors of 
the Pauli matrices a x , a y , and a z . Any unitary operation preserves the angles between the 
axes of the transformed operators, so we can redefine our coordinates to have a new set of 
Pauli matrices. We can then say that the structure of the MUBs remains invariant under 
any unitary transformation. This is akin to saying that only one MUB structure exists in 
the two-dimensional Hilbert space. 

The same result also holds for two qubits, although, in addition, the extra feature of 
entanglement appears. Several methods have been presented for the explicit construction 

of MUBs 0, lH H B H B Q- Here we wil1 onl Y focus on one of them > which is based 
on the use of the finite Fourier transform, employing the operators a x , and a z and tensor 



products 



(Because this work attempts to delineate the structure and interrelation 



between MUBs, and not their explicit mathematical construction, we omit such a discussion 
and direct interested readers to the aforementioned work.) If we follow the algorithm in [47], 
we get a table with five rows of three mutually commuting (tensor products of) operators, 
shown in Table HI (We have suppressed the tensor multip ication sign in all the tables.) The 
table reproduces Eqns. (3.30) and (3.32)-(3.335) in Ref. 4l|- 

By construction, the algorithm guarantees that the simultaneous eigenstates of the oper- 
ators in each row give a complete basis, and each basis is mutually unbiased to each other. 
The number on the left enumerates the bases, while the number on the right denotes how 
many subsystems the bases can be factorized into. 

It is easy to see that the three first bases are fully separable (the three operators in each 
of the first three rows commute for each of the two subsystems, separately), and that the 
last two bases are not separable. In fact, their simultaneous eigenstates are all maximally 
entangled states. We call this MUB construction with three (bi-) separable and two nonsep- 
arable MUBs a (3,2) construction. This is equivalent, under local unitary transformations, 
to the construction given in Fig. 1 of Ref. jsj ■ 
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1 




la z 


a z t 


2 


2 


&x&x 


ta x 


a x t 


2 


3 




lOy 


Oyl 


2 


4 


a x a y 


CTzCTx 


(Ty&z 


1 


5 


OyO X 


a z a y 


&X&Z 


1 



TABLE I: Five sets of three operators defining a (3,2) MUB. 

The algorithm imposes several characteristic features of the operator table. The table 
is composed of binary tensor products of the four operators a z , a x , a y = icr x a z = —ia z a x , 
and 1 = a 2 = a 2 = a 2 . In all, there exists 4 2 = 16 combinations of such products, but 
the operator i ® i must be excluded because it commutes with every operator in the set. 
Each of the 15 = 3-5 remaining operators is represented once in the table above. Moreover, 
the table is uniquely defined by the four entries in the two first columns of the first two 
rows. All other operators O rjC are determined by the relations O rc = O riC _ 2 O riC _i, and 
O r ^c = 2)C Oi )C+r _3 for r > 2, where the indices r and c denote the row and the column of 
the operator, respectively, and must be taken modulo four. 

Noting that each separable basis (i.e., the first three rows) has two eigenoperators con- 
taining the identity, that a nonseparable basis cannot have any eigenoperator containing the 
identity, and that there must be six entries, 1 Cg) a x , 1 ® a y , . . . , a y ® 1, a z (g) 1, containing 
the identity in the table, we can conclude that the (3,2) is the only possible construction in 
this space. That is, any unitary transformation, local or nonlocal, will yield an isomorphic 
table with respect to the separability, except, perhaps, for some row permutations. 



III. MUB STRUCTURES FOR THREE QUBITS 



Lawrence, Brukner, and Zeilinger [38] have shown explicitly two different sets of MUBs 
in the case d = 8. One of them has three fully separable bases (every eigenvector of these 
three bases is a tens or p roduct of states embedded in the Hilbert space of each single qubit) 
and six GHZ bases [48]. The other structure has nine sets of bases with eigenvectors where 
one qubit can be factorized and the ot 
we follow again the algorithm in Ref. 

Table ILTl is equivalent to the first MUB construction demonstrated in this space by Fields 



ler two qubits are in a maximally entangled state. If 
3 we get Table HH 
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a z tt 


Ha z 


la z t 


& z ta z 


la z a z 


&Z&Z&Z 


a z a z l 


3 


2 


a x ti 


ia x i 


iia x 


a x a x l 


\a x a x 


x @ x x 


a x TLa x 


3 


3 


a y it 


ta x a z 


la z a x 


&y&x&z 


\a y a y 




Gy&z&x 


2 


4 


a x la z 


Ifryt 


a z la y 


O'xd'yd'z 


& z a y a y 


&y v y x 


a y ta x 


2 


5 


a x a z t 




la z a y 


OyOyGz 


Vz&y&x 


GyV% @ x 


a x ta y 


1 


6 


a y ta z 


la y a z 


a z a z a y 


O-yfryl 


<?z<?x<?x 


@x@x@y 


O'xO'zO'x 


1 


7 




a z a y a z 


O'z^zO'x 


a y a x l 


la x a y 


@x®y@x 


<T y \a y 


1 


8 


dy& z d z 


& z a y t 


a z la x 


O'xO'x^z 


ta y & x 


@y® x® y 


O'x^zO'y 


1 


9 


a y a z l 


a z a x l 




a x a y l 


(J Z X 0y 


@x®y® y 


O y G z G y 


2 



TABLE II: Nine sets of operators denning a (2,3,4) MUB. 

and Wootters or to the example 2 of Sec. 5 in Ref. Q|, in that it has two fully (that 
is, tri-) separable bases (marked with a 3 in the rightmost column), three biseparable bases 
(marked with a 2), and four nonseparable bases (marked with a 1). We will denote such a 
set of MUBs SIS db (2,3,4) structure. 

Lawrence, Brukner, and Zeilinger j3| have pointed out that two other constructions are 
possible, namely a (3,0,6) set of MUBs where three of the bases are fully separable, and the 
remaining six bases are nonseparable, and a (0,9,0) in which all the bases are biseparable. 
The corresponding operators are given in Figs. 2 and 4 of Ref. (^J. We would like to see 
how we can derive these bases from the ones in Table ITT1 To this end we use the controlled- Z 
operator 



V 



(3.1) 



^10 ^ 
10 
10 
-1 

This operator is unitary, nonseparable and, moreover, has the property that it is its own 
inverse, Z~ l = Z c , and its own conjugate. It commutes with the operators t ® a z , a z (g) 1, 
and o z ®b z . Let us first convert the (2,3,4) into an equivalent basis set by applying the local 
unitary (permutation) transformation x—>y-^z-^xto the leftmost qubit. The operator 



6 



performing this transformation (up to an overall phase-factor) is 

1 



( gi7r/ '4 gi37r/4 



gJTr/4 g-JTr/4 



(3.2) 



We also apply the permutation y «-> z to the middle and rightmost qubits. The corresponding 
operator is 

1 



1 i 
z 1 



(3-3) 



Applying the operator U P ®U C ®U C to Table ITT1 above, we are left with an equivalent operator 
table, Table HH still defining a (2,3,4) MUB. 
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ii&y 


idyi 


a x ta y 




X <Jy<7y 


a x a y i 
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a y ii 


ia x i 


iidx 


a y a x t 


la x a x 


VyVxOx 


a y la x 


3 


3 


a z ii 


ia x & y 


l&y(Tx 


O z d x Oy 


ia z a z 


a z a z a z 


ezey<J x 
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4 




ia z i 


a x ta z 


a y a z a y 


&x&z&z 




a z ta x 
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5 


d-y&yl 


@x@x@y 


ta y a z 


o z o z dy 


&X&Z&X 


&z&x&x 


a y la z 


1 


6 


a z ta y 


la z a y 


VxVy&z 


a z a z l 


&X&X&X 


OyO x O Z 


O-yO-y&x 


1 


7 


(T y <J y a y 


O x O Z Oy 


v x O'yU' x 


a z a x t 


la x a z 


d y o z o x 


a z ia z 


1 


8 


<J z <JyGy 


a x a z t 


&x^-&x 


&y&x&y 


la z a x 


&z&x&z 


a y a y a z 


1 


9 


a z a y l 


& x o x \ 


Ha z 


a y & z t 


&x&x&z 


<Jy<J Z (j Z 


o z &yd z 


2 



TABLE III: A local unitary rotation of Table 1 defining a "different", but isomorphic (2,3,4) MUB. 

Now we apply the nonlocal unitary operator 1 ® Z c to the three qubits (where, evidently, 
the operator is only nonlocal in the four- dimensional subsystem constituting the rightmost 
two qubits). The local transformation we performed above on the first qubit will of course 
not change the structure of the MUB, not even after subsequently applying the operator 
1 ®Z C . The reason we made this local transformation is only to facilitate a comparison with 
the construction in Ref. We note that Z c [i ® (T X )Z\ = Z c {i ® <J x )Z c = d z ® a x , and 

therefore Z c {a z ®b x )Z\ = i®a x - The corresponding transformations for the other product 
operators are 



a x ® 1 <-> 0x® &z, 
i® a y <-> a z ® a y , 



(3.4) 



a y ®t 



a y ® o z . 
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From these, the remaining relations 

fr x ® fr x «-> fr y ® fr y ^ 

fr X ® fry ^fry® fr X 

follow. Hence, applying this transformation to Table IIII| will result in Table IIVI From the 
unitarity of Z c it follows that all inner products between the eigenstates of the simultaneous 
eigenvectors of the operators in the same or in different rows of the two tables will be 
identical. We can therefore be confident that Table IIVI corresponds to a set of MUBs. 
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a x ii 


idzCTy 


tfryfrz 


PxCzPy 


\a x a x 


@x®x®x 


CxPyVz 
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frytl 


ta x a z 


ia z frx 


fryfrxfrz 


i-fryfry 


OyOyOy 


fryfrz&x 


2 
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a z tt 


tfryfrx 


la x a y 


Ozfryfrx 


tfrzfrz 


d z o z frz 


a z fr x fr y 
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fryfrzfry 


lfr z l 


frxtfrz 


fryifry 


frxfrzfrz 


a z la x 


frzfrzfrx 
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fryfryfrz 


frxfryfrx 


ia y i 


a z ta y 


a x ifrx 


d z fryfry 


frylfrz 
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a z a z a y 


iifry 


a x a y l 


fr z a z t 


®x®y®y 


a y a x t 


®y@x@y 


2 
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fryfrxfrx 


a x tfr y 


@x®x®y 


frzfrxfrz 


ia x t 


&yta x 


fr z tfrz 


2 


8 


frzfrxfrx 


frxfrz't 


Oxfrzfrx 


aydyCx 


ii&x 


frzfrx 1 


fryfry\ 
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a z a y a z 




iifr z 


fryfrzl 


a x (Jxl 


a y a z a z 


a z cr y l 


2 



TABLE IV: Nine sets of operators defining a (0,9,0) MUB. 



However, this set of MUBs represents a different entanglement structure, because here 
every basis is biseparable. In our nomenclature it is a (0,9,0) set. In fact it is the same table 
(with some rows interchanged) as Fig. 4 in Ref. 

We can now continue and apply the operator Z c <g> 1 to Table llllj. above. Again the set 
of simultaneous eigenstates of the operators in any row will define a complete basis, and the 
set of bases will form a MUB. The result can be seen in Table fVl 

This yields a (1,6,2) MUB. That is, only one of the bases is fully separable. This 
construction is novel: it is neither a structure of the Fields' and Wootter's type, nor is it 
one of Lawrence et al's two structures. We note from the tables above, that there are nine 
operators containing two identity operators and twenty seven containing a single identity 
operator. In each of s sets of operators defining a fully separable basis (i.e., in each of s 
rows), there are three entries with two identity operators each. In each such set there are also 
three entries with a single identity operator. Each of b operator sets of defining a biseparable 
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y y 


~y ® x ® V 

y *" y 


&y& X t 
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a z a x t 


iia x 
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(Tz&x&x 


O x Oy(J x 
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b z CT y d x 


1a x a y 


ta z a z 


d'zd'zd'z 


\Oy(J X 


2 
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ia z i 




Oy\dy 


& x ta z 


d'z&zd'x 


& z ta x 
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d'z&z&y 


a x la x 


\& x (j x 


b y cr z a z 
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a z ta y 


la z ay 


VyVx&z 


& z fr z t 


@y@y@x 


X 0y& Z 


&X&X&X 
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@x®x®y 


a x ta y 


@y@x@x 


i& x i 


O'zO'x^z 


& y ta x 


a z ia z 
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a x it 


Vx<?z<?x 


®x ®y ®y 


ta z a x 


la x a z 


Vx&xVz 
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ia y i 


(Tydyl 


ii<5- 2 




bydyd z 


& y ta z 


tb y a z 


3 



TABLE V: Nine sets of operators denning a (1,6,2) MUB. 

basis contains one operator with two identities, and three operators with a single identity. 
Finally, the n sets of operators defining nonseparable bases contain no operators with two 
identities, and three operators with a single identity. In consequence, we have the equations 



3s + 6 = 9, 



3(s + b + n) = 27, 



(3.6) 



for all nonnegative integers (smaller or equal to nine), which yield the four solutions 
{(2, 3, 4), (0, 9, 0), (1, 6, 2), (3, 0, 6)}. We conclude that, so far, we have derived explicit con- 
structions for the first three structures, and have one more left to construct. 

Before doing that, we make a small digression and note that for three qutrits 
(4^]], similar considerations lead to the conclusion that in this 2 7- dimensional 
Hilbert space, with 28 MUBs, there exists five MUB structures, namely the set 
{(0, 12, 16), (1, 9, 18), (2, 6, 20), (3, 3, 22), (4, 0, 24)}. 

Now we return to the three qubit space. The last possible (3, 0, 6) structure can be built 
up in the following way: Take Table IVl and perform the transformation i/^zon the leftmost 
two qubits and the transformation x <-> y on the rightmost qubit. As we are only going to 
transform the two rightmost qubits in a nonlocal fashion, the transformation on the leftmost 



qubit is only to yield a table that is identical to one of the constructions in Ref. 



33. The 
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operator performing the transformation x «-> y is 

^ / e^/ 4 \ 
^ e-^/ 4 J ' 

The result of the local transformations £/ c £g> £7 C ® U r is shown in Table IVII 
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la z t 


a z a z t 


ita z 
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a z a z a z 


a z la z 


ta z a z 


3 



TABLE VI: A table, isomorphic to Table 4, defining a (1,6,2) MUB. 

If we subsequently transform Table IVII with il ® Z c , we get the last of the possible 
MUB constructions in this eight-dimensional space, which is reproduced in Table IVTT1 This 
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TABLE VII: Nine sets of operators defining a (3,0,6) MUB. 
table is a (3, 0, 6) MUB. It is, in fact, exactly the same (with some rows permuted) as in 
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(3.7) 



Fig. 2 of Ref . |38[ . The possibilities in the eight- dimensional space are now exhausted: No 
other MUBs with different entanglement structures can be constructed. We summarize the 
operational relationship between the different MUB structures one can construct in Fig. ^ 



Table II 

(2,3,4) 



Up®U c <g>U c 



Table III 

(2,3,4) 



1 ®Z C 



Table IV 

(0,9,0) 



Z c 8 1 



1 ® z c • U c S> U c g> u r 



Table V 
(1,6,2) 



Table VII 

(3,0,6) 



FIG. 1: The operational relationship between the different MUB constructions 



IV. MUB STRUCTURES FOR FOUR QUBITS 

With four qubits, the MUBs can take one of five different forms with respect to their 
separability. We have fully separable bases, triseparable bases (2x2x4), two kinds of 
biseparable bases (one that factors 2x8 and the other that factors 4x4), and finally 
nonseparable bases. If we follow once more the algorithm given in Ref. [47], but write 
explicitly only the first four columns of each basis to save space, we get the following table: 

The remaining eleven columns of the table can be generated through the relation T)C = 
O r)C -40 ryC -i. The 2 in the last column indicates a basis biseparable in a 2 x 8 space, while 
2B indicate a basis biseparable in a 4 x 4 space. This basis will be denoted as a (2,0,4,2,9) 
MUB, referring how many of the bases that are fully separable, triseparable, biseparable (in 
a 2 x 8 and in a 4 x 4 space, respectively), and nonseparable. We now apply the operator 



which first locally rotates the two rightmost qubits so that a y <->• d z and then entangles (or 
disentangles) the same two qubits. The result is Table llXl 

Before venturing further, it is instructive to see how many different MUB structures there 
are in the four-qubit space. Again we use the fact that identity operators play a special role 
in defining the separability of the bases. Table |X| shows how many products of three, two, 



l(g)l(g)Z c -t(g)l(g)U c (g)U ( 



(4.1) 
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TABLE VIII: Seventeen sets of four operators denning a (2,0,4,2,9) MUB. 

and single identity operators define a basis of a certain kind. The ordering of the multiple 
identities plays no role. 

Solving the three equations for the different number of identity operators, we find six- 
teen different MUB structures: (3,0,0,2,12), (2,0,4,2,9), (2,1,2,2,10), (2,2,0,2,11), (1,0,8,2,6), 
(1,1,6,2,7), (1,2,4,2,8), (1,3,2,2,9), (1,4,0,2,10), (0,0,12,2,3), (0,1,10,2,4), (0,2,8,2,5), 
(0,3,6,2,6), (0,4,4,2,7), (0,5,2,2,8), and (0,6,0,2,9). Of these, we have explicitly given the 
tables for the (2, 0, 4, 2, 9) and the (0, 4, 4, 2, 7) structures. Deriving the transformations be- 
tween any two of the sixteen structures goes beyond the scope of this paper. However, note 
that applying the operator 1 ® Z c <g> i t o the entries of Table I1X.I will yield a table with the 
MUB structure (0,3,6,2,6). Moreover, using instead the operator 

Z C ®\®\-U C ®U C ®\®\ (4.2) 
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TABLE IX: Seventeen sets of four operators defining a (0,4,4,2,7) MUB. 

in each entry of Table HXl will yield a (1,2,4,2,8) MUB structure. Since we know that a 
sequence of C-NOT or controlled- Z operations, together with local unitary rotations, suffice 
to make any entanglement transformation on qubits, it is clear that similar transformations 
will yield the whole set of different MUB structures, starting from the Table IVlTTl 

V. CONCLUSIONS 

For one and two qubits, all possible MUBs have the same structure with the respect 
of entanglement. For more qubits, the situation is more involved: four different MUB 
structures appear for three qubits and seventeen MUB structures for four qubits. The 
difference between these structures lies in how the bases are entangled. For both three and 
four qubits, MUBS exist that have 3, 2, 1 and no fully separable basis set(s). For three 
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TABLE X: The separability of the bases (left column) and the number of triplets, pairs, and single 
identity operators contained in the basis defining operator set. 

qubits it is possible to find one MUB that have no fully nonseparable bases. This is no 
longer possible in the four qubit case. 

In a quantum protocol relying on MUBs, the entanglement structure of the MUB is 
usually inconsequential. What counts is usually only the mutual unbiasedness, not the 
separability of the bases. Experimentally, however, it may be easier to generate one set of 
bases rather than another. Some of the bases can be generated locally, accessing each qubit 
separately. However, as we have shown, when several qubits are involved, most bases are 
entangled in one way or another, requiring joint operations on the qubits. In this paper we 
have tried to delineate the possible MUB structures for up to four qubits. The method we 
have employed can of course be extended to any number of qubits, although the complexity 
and variety of bases grows very rapidly with the number of qubits. 
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